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Recent numerical simulations of coalescing binary neutron 
stars conducted by Wilson, Mathews and Marronetti (WMM) 
show a rising central energy density of the stars as the orbital 
separation shrinks, i.e. the stars are individually crushed as 
they near coalescence. They claim this "star-crushing" effect 
is partially due to a non-linear, first post-Newtonian order 
enhancement of the self-gravity of each star caused by the 
presence of the other star. We present a concrete calculation 
which shows, within general relativity, first post-Newtonian 
order interactions with the other star leave the central energy 
density unchanged as the orbital radius shrinks. The results 
presented here are in sharp disagreement with the WMM 
claim. However, alternative gravitational theories, such as 
Brans-Dicke theory, can exhibit a small crushing effect in the 
binary constituents as they near coalescence. We show that 
the absence of the star-crushing effect at first post-Newtonian 
order is related to adherence to the strong equivalence prin- 
ciple. 

PACS Numbers: 97.80.Fk, 04.25.Dm, 04.40.Dg, 97.60. Jd 



Binary neutron stars spiraling toward their final coa- 
lescence under the dissipative influence of gravitational 
radiation reaction forces are the primary sources grav- 
itational waves which might be detected by interfero- 
metric gravitational wave detectors such as LIGO and 
VIRGO 10,0. Extracting the gravitational waves from 
the detector noise and using the information encoded in 
the signals requires a precise knowledge of the expected 
waveforms M . The primary method for studying the be- 
havior of these highly relativistic binaries has been the 
"post-Newtonian" approximation carried to high rela- 
tivistic order &I3 . Recently, Wilson, Mathews and Mar- 
ronetti (WMM) [7| have used numerical hydrodynamic 
simulations as an alternative method for studying the 
late stages of binary neutron-star inspiral. These simu- 
lations, based on the simplifying assumption that space- 
time is spatially conformally-flat (SCF) ||, obtain a se- 
quence of equilibrium configurations of the neutron stars 
at closer and closer orbital separation. Their numeri- 
cal results show that the neutron stars are individually 
compressed and the central densities begin to rise as the 
orbital separation is decreased. [See Mathews and Wil- 
son Figure 2.] Such a crushing-effect has never been 
reported in post-Newtonian calculations. This discrep- 
ancy forces a re-examination of the validity of both the 
post-Newtonian calculations and the numerical simula- 
tions. (See Eardley || for a discussion.) 



It is tempting to simply dismiss this unusual star- 
crushing effect as an artifact of the SCF constraint used 
in the WMM numerical simulations. However, WMM 
further claim that the star-crushing effect is partially 
due to a non-linear first post-Newtonian order (1PN) en- 
hancement of the gravitational potential. (See Mathews 
and Wilson Appendix A: "... the [first] post Newto- 
nian expansion does indeed exhibit an enhancement to 
the gravitational potential ...".) However, the 1PN met- 
ric is naturally spatially conformally flat; therefore such 
an effect would be independent of the SCF assumption. 
[See Eq.(l).] It is also tempting to dismiss the WMM 
claim on the grounds that Newtonian tidal-effects will 
stretch the individual stars, and thereby reduce the cen- 
tral densities of stars [Q. However, WMM have made it 
clear that their proposed crushing mechanism will swamp 
- by a factor of 10 4 - the rather small reduction in cen- 
tral density due to tidal distortion. [See Mathews and 
Wilson § Eq. (20).] 

There is actually a substantial body of literature 
on post-Newtonian effects for extended bodies in bina- 
ries that seemingly would have discovered a 1PN star- 
crushing effect if it were real. (See e.g. |fl],[l2).) There are 
also published arguments that show that enhancements 
of the local gravity due the presence of other gravitating 
bodies involve a 1PN violation of the strong equivalence 
principle. (See Will |§ Chapter 6.3.) However, WMM 
specifically claim a 1PN origin of the star-crushing effect; 
therefore we offer a specific rebuttal to the claim. 

In this paper we present a brutally straightforward, 
first post- Newtonian order calculation of the central en- 
ergy density of a compact star in a close binary. We 
solve the 1PN hydrodynamic equation (including non- 
linear gravitational interactions with the other star) for 
the density profile of the star. At intermediate steps in 
the calculation, relativistic corrections appear which give 
the illusion that each star is being compressed by the 
presence of the other star. In the end, for a first post- 
Newtonian expansion of general relativity, all the rela- 
tivistic corrections involving the other star cancel, and 
the true physical quantities (the central energy density, 
and the proper radius) are unaffected by the presence 
of the other star. Thus WMM's claim of a 1PN ori- 
gin of the star-crushing effect is false. The confusion 
can be summed up as follows: the potential [in this case 
(.9oo + 1)] is made "deeper" by non-linear relativistic fac- 
tors involving the other star; thus, when solving for the 
density profile, it appears that the star is being com- 
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pressed by the presence of the other star. However, the 
relativistic factors also re-enter the calculation as 1PN 
corrections to the spatial part of the metric. The spatial 
metric components alter the spatial line element in the 
calculation of the proper radius of the star, and they al- 
ter the proper volume element in the calculation of the 
conserved baryon number. In general relativity, at first 
post-Newtonian order, these effects exactly cancel, and 
the physical quantities are left unaffected by the other 
star. 

We consider two perfect fluid blobs (stars) with baryon 
masses M bA and MbB ■ For simplicity we assume that the 
stars are in a nearly circular orbit with coordinate radius 
R Q . We boost to a comoving, post-Newtonian coordinate 
system where the fluid in star-A is momentarily at rest, 
and the origin instantaneously coincides with the center 
of star-A. [Star-A is still being accelerated toward star B.] 
In this moving frame the parameterized, post-Newtonian 
(PPN) metric coefficients arc |p|,p^| 

5 oo = -1 + 2C/(x, t) - 2/3C/(x, t) 2 + 2( 7 + l)$ x (x, t) 
+2(3 7 - 2/3 -(- l)$a(x, t) + 2$ 3 (x, t) + 6 7 $ 4 (x, t) (la) 
gij = [l + 2 7 f/(x,t)]5 ij (lb) 
g oj = -(l/2)(4 7 + 3)VS-(x,<) - (l/2)W-(x,*) , (lc) 

where (3 and 7 are the PPN coefficients which charac- 
terize alternative theories of gravity. Solar system tests 
have constrained their values to be very close to unity. 
In general relativity fj = 7 = 1. Also notice the metric is 
spatially conformally flat (SCF) [gj. The sources for the 
potentials (the baryon density p, pressure P, and internal 
thermal energy pll) are measured in the comoving frame 
of the fluid. 

We now define and evaluate the potentials appearing 
in Eq.(l). (Details can be found in The Newtonian 
potential is 
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where L denotes the - as yet unknown - coordinate ra- 
dius of the star, R Q is the vector separation of the body 
centers, R = |R | is the constant coordinate orbital 
separation, x is the field point, r = |x| is the radial co- 
ordinate measured from the center of star-A. The sec- 
ond expression in Eq. (|2]) is valid only within the quasi- 
stationary star-A. We have also neglected tidal terms. 
The (R • x)-term gives rise to the bulk acceleration of 
star-A toward star-B. Over the extent of star-A, it has 
zero average radial (r) component; therefore it cannot 
contribute to the compression. We will neglect this term 
and other terms which clearly cannot contribute to the 
radial compression. The Newtonian potential can now 
be written 



U(r <L) = 
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Ro 



U A (r) 
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where U A (r) is defined as the "47r"-term in Eq. (|2|). For 
future reference, we compute the radial derivative 



U, r = — 
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p{r') r' dr' 



M bA (r) 



(4) 



The potential <i>i is the post-Newtonian contribution 
from the kinetic energy. As we have boosted to a frame 
where star-A is a rest, there is no bulk motion of the fluid 
in star-A to contribute to this term. We also neglect in- 
ternal fluid motion in star-A, as this actually decreases 
the central density (Hi. Thus we have 



$i(x,t) 



v 2 p(x',t)d 3 x' 



M bB 
R 
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where vb is the velocity of star-B measured in the frame 
where star-A is momentarily at rest. The potential <I>2 is 
the post-Newtonian contribution from the gravitational 
energy 



$ 2 (x,t) = J 



f U(x',t)p{x',t)d 3 x' 



(6) 



Proceeding as with previous terms we have 
M bB M A (r) 4tt f r 



$ 2 
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p A (r')U A (r') r' z dr' . 
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The origin of the WMM star-crushing claim is that the 
first term in Eq.(^) will conspire with the Newtonian 
gravitational pull in Eq.(Q) to enhance the compression 
of star-A. Noting that the WMM-factor (W 2 - 1) ~ 
MbB/Ro, we see that radial enhancement to the Newto- 
nian compressional force appearing in Eq. (0) is the force 
related to the potential energy shown in Mathews and 
Wilson @ Eq.(19) [AE GR , w 2(W 2 - \)GM 2 jL\. Thus, 
we are unambiguously addressing the term that WMM 
claim is the origin of the crushing-effect. 

The potentials $3 and $4 are similar to $1 with the ki- 
netic energy replaced by pll and P respectively. Finally, 
the potential can be evaluated using formulae in || The 
results: Vj = MbBV J B /R , and Wj — for circular orbits. 

We now look at at the spatial components of the Euler 
equation in our instantaneously comoving- moving frame 
of star-A. In this frame, the fluid is momentarily station- 
ary (u k = P t = for star-A), and we can write 



dP 

dr 



(p + pU + P) 



dr 



log V-.9oo 
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We have neglected a number or second post-Newtonian 
terms, as well as terms which do not contribute to the 
radial compression. 

Expanding the logarithm, substituting the potentials, 
and discarding second post-Newtonian terms, we have 
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dP 

dr 



pU, r = ep 



2(1 -p)UU, r +(3 7 - 2p+l)$ 2 , r 



+ $ 3 ,r +3 7 $4,r + [II + (P/ P )]U„ 



(9) 



On the left we have the familiar pressure gradient balanc- 
ing the Newtonian gravitational attraction. On the right 
we have all the first post-Newtonian corrections. We have 
inserted a dummy factor of e(= 1) to keep track of post- 
Newtonian terms. Multiplying by r 2 /p, differentiating, 
and using the expressions for the radial derivatives of the 
potentials, we have 



l-e(37-4^ + 3)(M 6B / J R ) 1 d 



r 2 dP A 



p dr 



PA 



4-7T r 2 dr 

?—lM A (r) 2 + PA [(37 - 4/3 + 3)U A + ^a] 

-(37- 



^ A dr^\ A pa ) 47rr 2 } 
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We have discarded terms of second post-Newtonian or- 
der, i.e. 0(e 2 ). On the left of Eq.(pX|) we have Newto- 
nian terms, as well as the post-Newtonian terms related 
to the presence of star-B. On the right we have only post- 
Newtonian self-interaction terms. The terms on the right 
will affect the matter configuration of star-A, but they 
would do so whether or not star-B is nearby. Our focus 
is on deriving the 1PN effect star-B has on the matter 
configuration in star-A. This can be done by setting the 
right hand side of Eq. (|l^) to zero and solving the homo- 
geneous equation. In a full solution to Eq.(10) any non- 
linear interplay between the post-Newtonian corrections 
on the left and right sides of Eq.(|lO|) would be 0[e 2 ], and 
could be neglected. Throughout the remainder of this 
calculation, we ruthlessly discard post-Newtonian self in- 
teraction terms. (A 1PN solution of Eq.([To|) - including 
the self interactions on the right side - will be presented 
elsewhere fl4| .) 

Substituting the polytropic equation of state 



P = Kp 1 



-l/r. 



(11) 



into the homogeneous equation, and defining the quanti- 
ties 



P= Pc 



(n + l)K 
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A= ap c 2n 
£= r/A . 



l-|^(37-4^ + 3) 



(12a) 
(12b) 

(12c) 
(12d) 



we obtain the dimensionless Lane-Emden Jig] equation 
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[We restrict our attention to non-singular (n < 3) solu- 
tions fl5| .1 The solutions 9(£;n) of this non- linear differ- 
ential equation must satisfy the boundary conditions 



0(0; n) 



1 



?'(0;n) = 
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The first condition defines p c as the central baryon den- 
sity. The second insures there is no pressure gradient 
at the origin. The density profile of star-A can then be 
written 



p A (r) = p c e(r/X;n) n 



(15) 



Notice that pa{t) has an apparent 1PN dependence 
on the potential produced by the other star (M(,s/i2 ) 
through Eq.(|l2]c). The coordinate radius L of the stellar 
surface is the value of r where the pressure (and thereby 
the density) vanishes; it is given by 



L — A fi„ 



(16) 



where fi„ is the smallest positive root satisfying 
0(fm;n) = O. 

The total baryon mass can be obtained by integrating 



Mf,A= / p <i(proper volume) = / pu°y/—gd 3 x 
J A J A 



[ p[l + 3ejU]d 3 

J A 



I A 

Aira 3 f(n)p~ 



1+^(49-7-3) 



2 Ro 

= M bA [K, n, pc, (4/3 - 7 - 3)(M hB /R )} 

3-n 

= A-na f(n)p c 2n (in gen. rel. = 7=1 
= MbA[K, n, p c ] (in general relativity ), 



(17a) 
(17b) 

(17c) 

(17d) 

(17e) 
(17f) 



where f(n) = fi n |#'(Cin, n)\ is solely a function of the 
polytropic index n. We have discarded self-interaction 
terms and 2PN terms. (We have also set the dummy 
e = 1.) Notice, - and this is the point of the paper - 
in general relativity at 1PN order for a star with fixed 
baryon mass (M^a) cind a fixed equation of state (K and 
n), there will be no change in the central density as the 
binary orbital separation changes. The functional depen- 
dence on Mbs/Ro cancels out of Eq. ( |l7|) for general rel- 
ativity (or any gravitational theory with 4/3 — 7 — 3 = 0). 
The inclusion of the self interactions that we have been 
discarding would complicate the functional form shown 
Eq. fll7jc), but it would not affect the the functional de- 
pendence depicted in Eq. (|l7jf); therefore our no-crushing 
conclusion would hold even if we were to include the self 
interactions on the right of Eq ( |Io[ ) . 

The proper radius, Aa, of the star can be obtained by 
integrating 



A.= 



10 

Ln 



1 



n Mi 
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(4/3 - 7 - 3) 



3 — n R a 

= Ln (in general relativity /3 = 7 = 1) 



(18a) 

(18b) 
(18c) 
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where 



f MbA 
\47ra 3 /(n) 



(19) 



is the radius of a the isolated Newtonian star with baryon 
mass MbA- Notice, in general relativity at 1PN order, the 
proper radius of the star in a close binary is unaffected 
by the presence of the other star. 

In the Brans-Dicke theory of gravity (3 = 1 and I7— 1| 
0.001 0]; therefore the compression effect illustrated in 
Eq. (|l7]) and (jl^) will be small. However, even with these 
constrained values of (3 and 7, this effect is potentially 
larger than those due to tidal distortion. [See Mathews 
and Wilson § Eq. (20).] 

The 1PN crushing effect we have been discussing can 
be interpreted as a change in the local value of Newton's 
gravitational constant due to the presence of the other 
star in the binary. That is, Eq. (ft7jc) and Eq. (|Tg|b) can 
be obtained from the standard Newtonian formula 
by the replacing G = 1 with 



Gi ocal = 1 - (4/? - 7 - 3)(M bB /R ) 



(20) 



Such effects arise in gravitational theories that vio- 
late the strong equivalence principle. [See reference ||] 
Eq. (6.75).] 

Notice, we have only shown that there is no crushing 
effect at 1PN order in general relativity. As the con- 
ventional (pefect fluid) post-Newtonian expansion of the 
metric is naturally spatially conformally flat at 1PN or- 
der, we cannot conclude with certainty that the crushing 
effect that WMM see is an artifact of the SCF constraint 
used in the numerical simulations. Indeed, it is possi- 
ble that the crushing effect is real, and WMM's only sin 
was to incorrectly assign a 1PN origin to the effect. The 
question remains whether the SCF boundary condition is 
somehow responsible for (or amplifies) the crushing effect 
in the simulations p6| . 

The post-Newtonian expansion of the metric is not 
spatially conformally flat at second post-Newtonian or- 
der Jll]|l8| ; therefore a thorough 2PN investigation may 
clarify the difference between the conventional post- 
Newtonian analysis of the inspiral problem and the SCF 
numerical analysis [Q. WMM claim that the SCF con- 
traint is justified, because it removes only the radiation 
from the solution. This is not true. Radiation only en- 
ters the metric at 2.5PN order; thus altering the met- 
ric at 2PN order by imposing the SCF constraint will 
change the longitudinal fields that govern the periodic 
orbital motion and the individual body configurations. 
Although it is true that the WMM equilibrium solutions 
are exact solutions (up to numerical accuracy) of the Ein- 
stein field equations, they are solutions which differ from 
the actual inspiral metric at the 2PN level. In this formal 
sense, the SCF simulations are only as accurate as a 1PN 
simulation. 
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